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We derive general expressions at the one-loop level for the coefficients of the covariant structure of 
the neutrino self-energy in the presence of a constant magnetic field. The neutrino energy spectrum 
and index of refraction are obtained for neutral and charged media in the strong-field limit (Mw 2£> 
s/B ^> m e , T, fi, |p|) using the lowest Landau level approximation. The results found within the 
lowest Landau level approximation are numerically validated, summing in all Landau levels, for 
strong B ^> T 2 and weakly-strong B > T 2 fields. The neutrino energy in leading order of the 
Fermi coupling constant is expressed as the sum of three terms: a kinetic-energy term, a term 
of interaction between the magnetic field and an induced neutrino magnetic moment, and a rest- 
energy term. The leading radiative correction to the kinetic-energy term depends linearly on the 
magnetic field strength and is independent of the chemical potential. The other two terms are 
only present in a charged medium. For strong and weakly-strong fields, it is found that the field- 
dependent correction to the neutrino energy in a neutral medium is much larger than the thermal 
one. Possible applications to cosmology and astrophysics are considered. 
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I. INTRODUCTION 



There are many astrophysical systems on which the physics of neutrinos in a magnetic field plays an important 
role. Let us recall that proto-neutron stars typically possess very strong magnetic fields. Large magnetic fields 
B = 10 12 — 10 14 G have been associated with the surface of supernovas [| and neutron stars Q, and fields perhaps 
as large as 10 16 G with magnetars Even larger fields could exist in the star's interior. It is presumed from the 
scalar virial theorem £| that the interior field in neutron stars could be as high as 10 18 G. A magnetic field as this 
(~ 10 18 G) in the interior of a compact star will be larger in two orders than the chemical potential characterizing its 
quark matter density. 

Unveiling the interconnection between the star magnetic field and its particle current flows could shed new light 
to the question of the star evolution. For example, it is well known that neutrinos drive supernova dynamics from 
beginning to end. Neutrino emission and interactions play a crucial role in core collapse supernovae [f|. Their 
eventual emission from the proto-neutron star contains nearly all the energy released in the star explosion. Neutrino 
luminosity, emissivity and the specific heat of the densest parts of the star are governed by charged and neutral 
current interactions involving matter at high densities and in the presence of strong magnetic fields. Thus, a total 
understanding of the star cooling mechanism in a strongly magnetized medium is crucial for astrophysics. 

On the other hand, the explanation of large-scale magnetic fields observed in a number of galaxies, and in clusters 
of galaxies [|| seems to require the existence of seed fields of primordial origin Q . According to several mechanisms 
P|, strong primordial fields ~ 10 24 G could be generated at the electroweak transition. Even larger fields have been 
associated with superconducting magnetic strings, which would generate fields ~ 10 30 G in their vicinity if created 
after inflation 0- 

Were primordial magnetic fields present in the early universe, they would have had non-trivial consequences for 
particle-physics cosmology. For instance, as it is well known, oscillations between neutrino flavors may change the 
relative abundance of neutrino species and may thereby affect primordial nucleosynthesis (for a recent review on 
neutrinos in cosmology see ^3)- Therefore, if a strong magnetic field (Mw ^> \f~B ^> m e , T, fi, |p|, with Mw and m e 
the W-boson and electron masses respectively) modifies the neutrino energy spectrum of different flavors in different 
ways, a primordial magnetic field can consequently influence the oscillation process in the primeval plasma 



'Presently on leave at: Department of Mathematics, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 
02139-4307 



2 



The propagation of neutrinos in magnetized media has been previously investigated by several authors |12|-|15|. 
Weak-field calculations were done for magnetized vacuum in Refs. [13, [13, and at T ^ and fi = (fi is the 
electric chemical potential) in Refs. |14|. In the /i = case, as long as B < T 2 , both the field- and temperature- 
dependent leading contributions to the neutrino energy resulted of 1/M^-order ^3"0I- I n the charged plasma 
[ijj, the field-dependent terms were much larger, ~ but they vanished in the (/i — > 0)-limit. The weak-field 

results of papers [12|-|l5| led to think that the magnetic- field effects could be significant in astrophysics, because 
of the field- and /x-dependent terms of order 1/Mj^, but were irrelevant (~ 1/M^-order) in the early universe due 
to its charge-symmetric character (/i = 0). Hence, it has been assumed that in cosmology the main contribution 
to neutrino energy was the purely thermal term of order T 4 /Myy [H| . However, as we will show below, a strong 
magnetic field (Mw 3> >/B S> Tn e ,T,fj,, |p|) gives rise to a new contribution to the neutrino energy that is linear in 
the field, independent of the chemical potential, and that is of the same order (l/Mj^) as the largest terms found in 
the weak-field char ged - medium case. This new result can turn magnetic-field effects relevant for cosmology. 

In recent papers ^3 , US , we investigated the effects of a strong magnetic field on neutrinos in magnetized vacuum 
(i.e. with T = and /j, = 0). There, to facilitate the calculations in the strong-field limit, we extended the Ritus' 
Ep-eigenfunction method of diagonalization of the Green functions of spin- 1/2 charged particle in electromagnetic 
field [H|-[2iJ, to the case of spin-I charged particles 0|- This formulation, which is particularly advantageous for 
strong-field calculations, provides an alternative method to the Schwinger approach to address QFT problems in 
electromagnetic backgrounds |2l|. The use of the Ritus' method resulted very convenient to study the neutrino-self- 
energy in magnetized media, since it allowed to diagonalize in momentum space both the electron and the W-boson 
Green's functions in the presence of a magnetic field. Ritus' formalism has also been recently applied to investigate 
non-perturbative QFT in electromagnetic backgrounds |22| . 

From the above considerations it is clear that strong magnetic fields can play a significant role in a variety of 
astrophysical systems, and possibly also in the early universe. For these applications, the analysis has to be carried 
out in the presence of a medium. Thus, in the present paper we extend the results obtained in papers anc l S3 to 
include finite temperature and density, performing a detailed study of the effects of a strong magnetic field on neutrino 
propagation in neutral and charged media, and discussing possible applications to astrophysics and cosmology. 

We stress that in calculating the neutrino self-energy in a magnetized medium, we should consider, as usual, its 
vacuum and statistical parts. In this case the vacuum part depends on the magnetic field, and for strong fields it can 
make important contributions even at high temperatures T 2 > eB. The reason is that the vacuum and statistical 
terms have different analytical behaviors, due to the lack of the statistic ultraviolet cutoff in the vacuum part. This 
fact gives rise to a field-dependent vacuum contribution (I/M^-order) which is larger than the thermal one (1/M^- 
order), into the self-energy. Therefore, as shown in this paper, a strong magnetic field can become more relevant than 
temperature for neutrino propagation in neutral media. 

The plan of the paper is as follows. In Sec. ^ we consider the radiative correction to the neutrino dispersion 
relation in the presence of a constant magnetic field. We introduce the general covariant structure of the neutrino 
self-energy in the presence of an external field, and find the dispersion relation as a function of the coefficients of 
each independent term in the covariant structure. The general form of the found dispersion relation goes beyond 
any given approximation and medium characteristic and serves as a guidance for particular applications. The general 
expressions for the coefficients of the self-energy structures in the presence of the magnetic field are found in the one- 
loop approximation in Sec. IIIII The leading behavior in the l/M^ expansion of those coefficients are then calculated 
in Sec. IIVI in the strong-field limit (i.e. in the lowest Landau level (LLL) approximation) for neutral and charged 
magnetized media. These results are then used in Sec. to find the corresponding neutrino dispersion relations and 
the indexes of refraction in neutral and charged strongly magnetized media. In Sec. IVII the LLL-approximation is 
numerically corroborated by summing in all Landau levels and finding the values of the coefficients for parameter 
ranges corresponding to strong B ^ T 2 and weakly strong B > T 2 fields. Possible applications to cosmology and 
astrophysics are discussed in Sec. I Vlll Finally, in Sec. IVllil we summarize the main outcomes of the paper and make 
some final remarks. 



II. NEUTRINO SELF-ENERGY GENERAL STRUCTURE IN A MAGNETIC FIELD 

The neutrino field equation of motion in a magnetized medium, including radiative corrections, is 



where the neutrino self-energy operator XXp) depends on the parameters characterizing the medium, as for instance, 
temperature, magnetic field, particle density, etc. 




(1) 
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The operator JD(p) is a Lorentz scalar that can be formed in the spinorial space taking the contractions between 
the characteristic vectors and tensors of the system with all the independent elements of the Dirac ring. Explicit 
chirality reduces it to 

£ = *5> T= v ^ ( 2 ) 

where L,R = s(l ± 75) are the chiral projector operators, and is a Lorentz vector that can be spanned as a 
superposition of four basic vectors that can be formed from the characteristic tensors of the problem. In a magnetized 
medium, besides the neutrino four- momentum p^ , we have to consider the magnetic- field strength tensor F^, to 
form the covariant structures 

5^(p, B) = + bj>^ + cp^F^j" + idp^uj". (3) 

In @ we introduced the notations F pL , = r^rF^, and F ^ v — ■^B\ £ ^vp\F pX ■ In the covariant representation the 

magnetic field can be expressed as — \e^ vp \u v F p , where u M is the vector four- velocity of the center of mass of 
the magnetized medium. The presence of the magnetic field, and hence of the dimensionless magnetic field tensor 

F^ v and its dual F^, allows the covariant separation in (J3J between longitudinal and transverse momentum terms 
that appears naturally in magnetic backgrounds 

j*U = p^F, PF pi ,Y, i>x = P^h (4) 

Coefficients a, b, c, and d are Lorentz scalars that depend on the parameters of the theory and the approximation 
used. Notice that finite temperature and/or density would not introduce any new term into Eq. ©. This formulation 
is indeed the most general one for a magnetized medium. In magnetized space, even if T = fx = 0, the vector 
is present, since the presence of a constant magnetic field fixes a special Lorentz frame: the rest frame (on which 
Ufi = (1,0,0,0)) where the magnetic field is defined (u fl F pl/ = 0). Thus, the inclusion of a medium does not break 
any additional symmetry. To express ^2(p, B) in terms of m m only means, therefore, to make a change in the selected 
basis vectors 1 . 

The neutrino energy spectrum in a magnetic background can be found from Eq. Q as the nontrivial solution of 



det 



t + J2(P,B)]=0. (5) 



Using the covariant structure J3J), the dispersion relation (J5J takes the form 



2 2 (1 + bf - d 2 2 

Po = P3 + T , ; .2 zPx- (6) 

1 + a) — c z 



One of the main goals of this paper is to find the coefficients a, b, c, and d, in the one-loop approximation, and 
hence, the dispersion relations for different systems in the presence of strong magnetic-field backgrounds. 



III. ONE-LOOP NEUTRINO SELF-ENERGY 



Let us consider the one-loop corrections to the neutrino self-energy in the presence of a constant magnetic field. To 
leading order in the Fermi coupling constant the main field-dependent contribution to the self-energy comes from the 
bubble diagram (Fig.I-a) with internal lines of virtual charged leptons and W-bosons, and from the tadpole diagram 
(Fig.I-b) with virtual loop of charged leptons. 

In a neutral medium only the bubble diagram contributes, while in a charged medium both diagrams should be 
considered. However, only the bubble distinguishes between neutrino flavors, since the flavor of the internal charged 
lepton is directly associated to the flavor of the propagating neutrino (i.e. the one appearing in the external legs). For 
the tadpole, because the internal boson is the neutral Z-boson, the charged leptons corresponding to different families 
are linked to the same neutrino flavor. Therefore, in regards to magnetic field effects on neutrino oscillations, the 
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bubble contribution is the essential one. Thus, because the ultimate goal of our study is to elucidate the possible effect 
of a strong magnetic field on neutrino flavor oscillations, henceforth we only consider the bubble diagram contribution, 
represented in Fig.l-a, into the neutrino self-energy. 
The bubble contribution to the one-loop self-energy is 

Z(x,y) = t ±-Ry v S(x,y)'fG F (x,y) M v L (7) 

where S(x, y) and Gf(x, y) v M are the Green's functions of the electron and W-boson, in the presence of the magnetic 
field, respectively. Since both virtual particles are electrically charged, the ma gnet ic field interacts with both of them 
producing the Landau quantization of the corresponding transverse momenta |17|. [l8|. 




Fig. I. One-loop contributions to the neutrino self-energy in a magnetized medium, (a) Bubble graph - The dashed 
line represents the test neutrino, the solid line a charged lepton of the same family of the test neutrino, and the 
wiggly line the W-boson. (b) Tadpole graph - The dashed line represents the test neutrino, the solid line a charged 
lepton of any species, and the zigzag line the Z-boson. 

The electron Green's function, diagonal in momentum space at arbitrary field stren gth , was obtained in Refs. |19| , 
pfij by Ritus, using what has become an alternative method to the Schwinger approach |2l| to deal with QFT problems 
on electromagnetic backgrounds. In this method the electron Green's function in configuration space is given by 

S{ ?^ = ij£l* E * (*) -=1T-E q (v) (8) 
J7 v (27r) 7-9 + m e 

where <f M = (<7oi 0, <?2) 53), E q = 7°-E]7°, and the magnetic field has been specialized in the rest frame along the 
iJ-direction (i.e. given in the Landau gauge as Ajf* = BxiS^)- 

The transformation functions E q (x) in JSJl play the role, in the presence of magnetic fields, of the usual Fourier 
functions e lqx in the free case and are given by 

E q {x) =Y J E q „(x)&{(j), (9) 

with spin matrix 

A(o-) = diag(8 cr i,S a _ 1 ,6 al ,5 (7 _ 1 ), a = ±1, (10) 

and eigenfunctions 

E qa (x) = N(n)expi(q-x)D K (p) (11) 

where N(k) — (Air \eB\)* j\fn\ is a normalization factor and D K (p) denotes the parabolic cylinder functions |23 with 
argument p — y/2 \eB\(xi — -^) and positive integer index 

K = K (l,a) = l + ~-±, « = 0,1,2,... (12) 

The integer I in Eq. (|12|l labels the Landau levels I — 0, 1, 2, .... 

The electron momentum eigenvalue in the presence of the magnetic field, q, is given by 



<7 M = (?o, 0, -sgn (eB) ^2 \eB\ I, q 3 ^ 



(13) 
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Ritus' technique for spin-1/2 particles was recently extended to the spin-1 particle case in Refs. [Tg- In this 
case the corresponding W-boson Green's function can be diagonalized in momentum space as 

G F (k, fc% v = (2ir) 4 ?W (fc - k') _ 2 ^ - (14) 

k 

where S^'(k — fc') = i5 3 (fc — k')5 nn ', and the W-boson momentum eigenvalue in the presence of the magnetic field, fc, 
is given by the relation 



kl + k\ + 2(n - l/2)eB, n = 0,1,2,... 



(15) 



with n denoting the W-boson Landau levels. Considering the mass-shell condition fc = — M^. in Eq. (|15|) . we can 
identify at eB > Myy the so called "zero-mode problem" l24l-|2fij . As known, at those values of the magnetic field a 
vacuum instability appears giving rise to W-condensation 25] . In our calculations we restrict the magnitude of the 
magnetic field to eB < Myy> so to avoid the presence of any tachyonic mode. 
The W-boson Green's function in configuration space is then 



k 



'GO 



where the transformation functions are 



with 



T a ktl (x) = P% [T k {x)V x p- 



pa 

7_ V2 



/a/2 \ 
110 
i — i 
V V2 J 



and 



7}=0,±1 

The matrices are explicitly given in terms of the W-boson spin projections rj by 

= diag(5 Vt0 , 5 Vt i,5 n ,-i,5 Vj o), V = 0, ±1 

and the eigenfunctions J- k r) are 

T kr ,(x) = N(K')expi(k-x)D K ,($ 



(16) 



(17) 



(18) 



(19) 



(20) 



(21) 



with N(k') a normalization factor similar to that appearing in i|ll[l . and D K >(£) denoting the parabolic cylinder 
functions with argument £ = y/2 \eB\ (xx — k^jeB) and positive integer index k! given in terms of the W-boson 
Landau numbers n and spin projections 77 as 



k' (n, rj) = n — r\ — 1, k = 0, 1, 2, . 



(22) 



As the neutrino is an electrically neutral particle, the transformation to momentum space of its self-energy can be 
carried out by the usual Fourier transform 

(2vr) 4 S (4) (p - j/)£(p, B) = J d 4 xd 4 ye- l ( p - x ~ p '- v ^(x, y) (23) 

Substituting with Q, © and JTJjJ in the RHS of Q we obtain 
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4 -t(p.x-p' .y) I R 



1" 



' E q {y) 



j.q + m 



r 



V— 



d 3 k r^A^L"^) 

(2tt) 4 fc 2 +M£, , 



(24) 



The coefScients of the covariant expression for E(p, B) can be found using the explicit form l|24(l in the covariant 
expression 



(2^) 4 5^ (p - p')S(p, B) = aV|| + &V± + cYF^Y + idYF^l" 



(25) 



Introducing in l|24|) the transformation functions 2? g and M , and taking into account the following properties of 
the spinor matrices 

A(±) t = A(±), A(±)A(±) = A(±), A(±)A( T ) = 
7"A(±) = A(±) 7 II ! 7 ± A(±) = A(T)7 ± ) 



LA(±) = A(±)Z, RA (±) = A (±) R 



(26) 



where the notation 7" = (7°,7 3 ) and 7 1 - = (7 1 ,7 2 ) was used, we find for the coefficients appearing in l|25|l the 
following general expressions in the one-loop approximation (Henceforth we consider sgn(eB) > 0, to simplify the 
notation.) 



l,n 



£ / / d 'y 



d 3 q f d 3 k 



-ifjp.w—p'.y) 



(2^) 4 J (2^) 4 (f + m 2) (t + M^) 



{(•20 +Q3)(P3 ~ Po)In-2,l-l(x)In-2,l-l(y) + (<?3 - Qo) (?3 + Po 0*0^n,Z (j/) } 



(27) 



d 3 9 /■ d 3 fc 



-i(p.x-p' .y) 



(27T) 4 7 (27T) 4 (^ +m 2)^ 2 +M ^ 



{(p 2 + ipi)q 2 I n -l,l (x)In-l,l-l(v) + (P2 - ^l)^-^-!^-!^)^-!,;^)} 



(28) 



V / j4 / j4 I f d ^ e -Kp-x-p'-y) 

2p~ 2 ^ 



l,n 



{2-Kf J {2vf (f +m 2)(t + M* W 



{((I0 + q3)(P3 ~ Po)In-2J-l{x)In-2J-l(y) ~ (?3 ~ ?o) (?3 + Po)if»,J (^)-C.i (#) } 



(29) 



2^1 



d 3 q /• d 3 fc 



-i(p.x-p' .y ) 



(2nfJ (2 7 r) 4 (^ + m 2)^ 2 + M 2 / 



{(^2 + »Pl)?2- r n,-X 1 i( a: )-C-l 1 J-l(j/) - (P2 ~ *Pl)52-fn-l,J-l(^)-C-l,j(j/)} 



(30) 
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where 

I n>l (x) = 2nv^exp\ix- (k + q)}(p n (Z/V2)ipi{p/V2) (31) 

and tp m [x) are the orthonormalized harmonic oscillator wave functions, defined in terms of the Hermite polynomials 
H m (x) as 

2 — rn 1 2 

<p m (x/V2) = H m (x/V2)exp-(x 2 /4) (32) 

Expressions l|27() - (|30|l . when substituted in l|25[) . give the general formula for the one-loop neutrino self-energy in a 
constant magnetic field of arbitrary strength. Note that, in this approach, the W-boson/magnetic- field interaction is 

2 

kept in the poles of the self-energy operator through the effective momentum k , as well as in the harmonic oscillator 
wave functions (^ rl (f/v / 2)- 



IV. NEUTRINO SELF-ENERGY IN THE STRONG-FIELD LIMIT 



Since in the strong-field limit (Mw ^ V B ^> m e ,T,fi, |p|) the gap between the electron Landau levels is larger 
than the rest of the parameters entering in the electron energy, it is consistent to use the LLL approximation in the 
electron Green's function, while in the W-boson Green's function, because My/ \/B, we must sum in all W-boson 
Landau levels. 

The neutrino self-energy in a magnetized vacuum was found within this approximation in Ref. [Tt^ . In what follows, 
we extend that result to neutral and charged media, introducing finite temperature and density effects. However, 
since the vacuum contribution is always present in the statistical calculations of the self-energy, we summarize below 
the results found in Ref. |l7j . 

In the strong field approximation (1 = 0) the covariant structure of the neutrino self-energy reduces to 

(2^) 4 5^ (p - p')Z (p, B) = + dtfPF^i' (33) 

That is, b' and d! are zero and 

_ <*+»). E j d , x f d i y f ^^.^/^^ - - (34) 



i^fJ (2^) 4 ( ? | + mc )p + M^) 



(35) 



After summing in n and integrating in x, y, and k we obtain 



l9MP3 + P0) -e-^° I A„ r^TlJ^ (36) 



2pjj J ^[ql + m e }[(q\\-P\\) 2 -eB + M^] 

r 2 



Integrating in qu and considering that 3> eB,m 2 ,pf,, we obtain for the neutrino self-energy in the strongly 



magnetized vacuum 

Eo(p, B) = ao3*|| + coP^F^l" (37) 

with coefficients given by 

°° = -c^W$k. e ~ A '*° (38) 

Notice the 1/M^-order of the leading contribution in the LLL approximation. Using the zero-temperature weak- 
field results of Ref. we can identify the scalar coefficients of the general structure 10 for that case as 



a = 6 -0, c = --do ~ fjZ\¥T72 ( 39 ) 



3~ 6g 2 eB 

If we compare the vacuum results at weak- field, Eq. I|39|l , with those at strong field, Eq. (|38|l , we can see that in 
both cases the non-zero coefficients have field-dependent leading contributions of 1/M^-order. However, as we will 
show in Sec.[V] they will play very different roles in the neutrino energy spectrum, since they are associated to different 
self-energy structures. While the strong-field results i|38[) produce a magnetic-field dependence in the neutrino energy 
spectrum which is linear in the Fermi coupling constant, the weak-field results (|39l) generate a smaller second-order 
contribution. 
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A. Neutral Medium in Strong Magnetic Field 

In a neutral medium the self-energy operator can be written as 

S(p,fl) = E (p,S) + S T (p,S), (40) 

where So is the vacuum contribution (given in the strong-field approximation by Eqs. (|37fl - l|38[) ). and St is the 
statistical part which depends on temperature and magnetic field (in a charged medium the statistical part can also 
depend, of course, on the chemical potential). Our goal now is to find Sy(p, B) in the parameter range M\y ^> 
\f~B 3> m e ,T, |p|. Here we can also assume that the electrons are confined to the LLL. Thus, performing the same 
calculations as in 13411 - 13tj|) we arrive to a self-energy operator S(p, B) with a structure identical to (|33|) . and coefficients 

oo 

given by an expression similar to (|3fc>|) . but with the Matsubara replacement, J dqo — > (27ri//3) (?4 — ; 

m— — oo 

m = 0, ±1, ±2, ...), corresponding to the discretization of the fourth component (74 of the fermion momentum, obtained 
after introducing the Wick rotation to Euclidean space (qo — > 194), (3 being the inverse of the temperature (3 = l/T. 

After summing in m and extracting the vacuum part Ij36|l ■ we obtain that the coefficients corresponding to the 
thermal contribution St(p, B) are given by 

a' (T) = — Cq(T) = + M e-A/^gW ^ _ pl) 

, f qsjpj + 4 - ej) + 2»p 4 £i . . > (g3-»P4)(P4"£|+£i)-2ip4£l , , .1 

d<?3 l£i[(^ + e |-e?) 2 + 4 £ ^]^ (£l) + -^WTi5F /s(£2) l (41) 



where 



= itW' /s(£2) = 1-W (42) 



are the fermion and boson distribution functions respectively, with effective energies 



e x = ^qj + ml, e 2 = y/(q 3 - p 3 ) 2 + - 2eB (43) 

To do the integral in (|41|l we take into account the approximation Myy ^> \f~B > T ^s> m e , \p\. Thus, to leading 
order we can neglect the contribution of the boson distribution, since it is damped by the exponential factor e~ Mw ^ T . 
Integrating in momentum and expanding in powers of 1/Myy, we obtain, up to leading order, 

Er(p, B) = a (T)^n + coiVF^Y (44) 

"°< T > - -« tr > = ^sf'" 1 '"" < 45 > 

Hence, the thermal contribution is one order smaller in powers of than the field-dependent vacuum contri- 

bution H38fl . Notice that in the strong-field approximation the thermal contribution Q45f) has the same 1/M^-order 
as the corresponding thermal contributions at weak 01 an( A zero field. 

We should call attention to the fact that when calculating the thermal contribution to the e-W bubble by using the 
Schwinger proper-time method pl| , some authors |14| have considered the following expansion 



aik - U] >" "^ + ^ A M 4 A ^ +0( ^ ) (46) 
Ivl w w w 



within the W-boson Green function in the presence of a constant magnetic field 



where 



Ou-.u ),,,.= <-(,'.,) / 0e*( l -»lG(P) M „, (47) 



{x,y) = expi^ yil F^x v (48) 
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is the well-known phase factor depending on the applied field |27|. and A(fc) is the energy-momentum transfer. 

Based on this expansion, it has been argued that the non-local lower contribution to the bubble diagram should 
be of order 1/M^, and that the local (l/M|L)-order term only contributes at finite density. While this argument is 
correct when applied to the statistical part of the self-energy operator (it explains the l/M^ order appearing in the 
thermal coefficients ao(T) and cq(T) in l|45|)l. it is not valid for the self-energy vacuum contribution. For the vacuum 
(T = 0, /i = 0) contribution the situation is different, since we cannot neglect the internal momenta as compared 
with Mw in the poles of the W-boson Green's function, due to the lack of the ultraviolet cutoff which is present in 
statistics. Thus, when calculating the vacuum part, the first term in (|46|1 should be replaced by 5 fJlU /(M 2 v + k 2 ). Such 
momentum dependence in the pole of the W-boson Green's function makes a nontrivial contribution of (l/M^)-order 
into the non-local part of the bubble diagram. Observe that the W-boson Green's function obtained in our approach, 
Eq. I|14p. only has such structure indeed, and that the vacuum results at strong l|38f) and weak (|39fl fields give rise to 
non-local (l/M^ / )-order contributions to neutrino self-energy. 



B. Charged Medium in Strong Magnetic Field 

We consider now a medium with a finite density of electric charge. As usual, a finite density is reflected through 
the shift qo — *• qo — A* in the momenta of the particles with non-zero charge associated to /i. The chemical potential, 
depending on the value of the charge density, plays the role of an external parameter. In stellar medium, for example, 
the electric charge is fixed by the net charge of the baryons, which, due to their large masses, are usually treated as 
a classical background. Here the following comment is in order, in electroweak matter at finite density there exists 
the possibility to induce additional "chemical potentials" ■ These chemical potentials are nothing but dynamically 
generated background fields given by the average (Wq ) and (Bq) of the zero components of the gauge fields. They 
are known to appear, for instance, in the presence of finite lepton/baryon density 28] (for recent applications related 
to this effect see Refs. 29] ). In our case, a possible consequence of the condensation of such average fields could be 
the modification of the effective chemical potential of the W-boson, which then might be different from the electron 
chemical potential. Nevertheless, that modification, even if present in the case under study, will have not any relevant 
consequence, as the W-boson contribution will be always suppressed by the factor e ~ M w/T _ 

In the strong-field approximation (i.e. after considering the LLL approximation in the neutrino self-energy l|25|) - (|30|) ') 
the non-zero coefficients of the neutrino self-energy in the charged medium are given in Euclidean space by 

Pp\ M) 2 +e?][(«4 -Pi) 2 + 4Y 

oo 

This expression is similar to Eq. (|36fl after taking into account the Matsubara replacement J dqo — > (i / (3) (<Z4 = 



m— — oo 



m ,a ' m = 0j ...), with the additional change q& — ► q%, where g| = q^ + ifi. 

After carrying out the temperature sum in (|49|l . and subtracting the vacuum part H36|) . we arrive to 

a {T, H) = -c (T, n) = —2 e S [ >{p-p) dq 3 

2pii J 



(g3 + gl) 

ei[(iei ~Pi) 2 +e\] 



(93 ~ £1) 

£l[(«£l +P4) 2 - 



(g 3 + e 2 - ip 4 ) (g 3 - e 2 - « P4 ) 

s 2 [{is2 + P4) 2 + s 2 } JBK 2) e 2 [(ie 2 - P4 ) 2 +£ 2 } JBy 2U V ' 

where 

= i + X^r f »^ = i-Xw> (51) 

are the fermion/anti-fermion and boson/anti-boson distributions respectively. Here, we obtain an equal energy split 
for both distributions since both charged particles have the same electric charge and thus are characterized by the 
same chemical potential. 
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Assuming the approximation Mw ^S> \ B ^> fx, T, m e , |p|, Eq. H50|) is reduced in leading order to 

a' (T,») = -c^(T, M ) = -^Vfe+^) [jV _ _ N + ]e -A/2eB 6( 4 ){p _ pl) (52) 

where the notation N^ 1 for the electron/positron number densities in the LLL 

= \eB\ J 0/^(£i). (53) 

was introduced. 

It is interesting to consider the situation where /i > T, which is a natural condition in many astrophysical environ- 
ments. In this case Eq. I|52|) becomes 

a' (T, M ) = -c' oM = Z^V^W+^ ^/a^^,^ (54) 

lvl wP\\ 

After analytic continuation to Minkowski space, we find that the leading contribution to the statistical part of the 
self-energy in a strongly magnetized charged medium is 

S T)/1 (p,S) = a (T,^n + c Q (T,fi)p»F^Y (55) 

ao(T, „) = -co(T, „) = 2(2 , )2M ^| « (56) 

The apparent dependence on the momentum components poi P3 of the coefficients ao(T, /z), cq(T,/i) is deceiving, 
as we can easily rewrite Eqs. (|55|) - (|56|l in a more convenient way using the base formed by the four- velocity and 
the covariant magnetic field vector = \e^ up \u v F pX , as 

ZtAP' B )=MT,^ + MT,h)$ (57) 

From H57fl - H58[) we see that the leading term in the statistical part of the self-energy in the charged medium in strong 
magnetic field is independent of the momentum and proportional to 1/M^, thus of larger order than the statistical 
contribution of the neutral case (Eq. 1)45(1 ^. 

Comparing our result (|57|l - l|58|l with those found at weak field ^5|, we see that the coefficient of the structure ^ 
has similar qualitative characteristics. That is, both c coefficients are linear in /i and have the same order, 
Nevertheless, the coefficients of the structure t/l significantly differ in the strong and weak cases. In the weak-field 
limit a only depends on the density, thus its leading contribution is proportional to /i 3 (this is the characteristic term 
of the MSW effect H^|). However, for fields larger than fi, the relation /iB ^ fi 3 holds, and the leading becomes 
proportional to \iB, as in Eq. (|58|l . Thus, in the strong field case it is the field the parameter that drives the main 
dependence of both structures in St,^. 

V. NEUTRINO INDEX OF REFRACTION IN STRONG MAGNETIC FIELD 

The neutrino index of refraction, defined by 

n =wh (59) 

where p and E are the neutrino momentum and energy respectively, is a quantity that plays a significant role in 
neutrino flavor oscillations in a medium 26]. A flavor-dependent index of refraction would enhance the periodical 
change between different flavors of neutrinos travelling through the medium. 
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To find the index of refraction we need to know the dispersion relation E = EQ p |). In the presence of a medium 
the energy spectra of the massless left-handed neutrinos are modified due to their weak interaction with the particles 
of the background. In a magnetized medium, even though the neutrinos are electrically neutral, they feel the magnetic 
field indirectly through their interaction with the charged particles of the medium whose properties are affected by 
the field. 

As discussed in Sec^J in a magnetized medium the general form of the neutrino dispersion relation is given by 
Eq. (JjjJ, with the self-energy depending on the medium characteristics. As follows, we present the neutrino energy 
spectrum and index of refraction for magnetized neutral and charged media. 

A. Neutrino Index of Refraction in Strongly Magnetized Neutral Medium 

In the strongly magnetized neutral medium the dispersion relation is obtained from 

det [i> + £ (P, B) + St(jj, B)} = 0, (60) 

with So(p, B) and Ht(j>,B) explicitly given in Eqs. (p?7| - ll3%|l and l|H| - (|^5|l respectively. The solution of lf6"U)l in 

leading order in the Fermi coupling constant is 



Pa = ± 



2/o(°)(p x b; 



± |p| [1 - a (0) sin" a] (61) 



In (|6ip. a is the angle between the direction of the neutrino momentum and that of the applied magnetic field, and 
the coefficient = ao + clq(T) is explicitly given by 



a 



8M 2 , L (2^)2 3M 2 , 



]exp(-pi/efl), (62) 



From l|62(l it is clear that in the strong-field limit the thermal correction eto(T) to the neutrino energy is much 
smaller than the vacuum correction ao- As we will discuss in Sec lVIIll this result can be significant for cosmological 
applications in case that a strong primordial magnetic field could exist in the primeval plasma. 

The neutrino (+)/anti-neutrino (— ) energies are given by E± = ±po respectively. Substituting Eq. (|61|l into i|59|) 
we obtain for the neutrino (antineutrino) index of refraction in a neutral medium 

sin 2 a (63) 

Eq. (|63|l implies that the index of refraction depends on the neutrino direction of motion. The order of the 

anisotropy is g 2 j^-- Maximum field effects take place for neutrinos propagating perpendicularly to the magnetic 

field. That is, the maximum departure of the neutrino phase velocity from the light velocity c, occurs at a = ir/2. 
Notice that in the neutral magnetized medium the magnetic field effect does not differentiate between neutrinos and 
antineutrinos. 

Here the following comment is in order. If we consider the weak-field results l|39|) obtained in Ref. ^3 f° r the 
neutrino self-energy, in the dispersion relation ijBJ) we obtain 

£± = |p|[l + ^sin 2 «] (64) 

In (|64(l . as the energy depends onc§, we can see that the weak field produces a negligible second order correction in 
term of the Fermi coupling constant expansion (i.e. a 1/M^-order effect). It can be corroborated that the inclusion 
of temperature in this approximation also produces a second order correction |14| . 

Thus, we conclude that the strong field produces an effect qualitatively larger, which is even more important than 
the thermal one. We call reader's attention to the fact that the field-dependent vacuum contribution to the neutrino 
energy l|61l) has the same order in the Fermi coupling constant as the ones found in a charged medium at zero pt| 
and weak |l5j fields. 
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B. Neutrino Index of Refraction in Strongly Magnetized Charged Medium 

Now we consider the effect of a strong magnetic field in neutrino propagation in a charged medium. In this case 
the dispersion relation is 

det [i> + S (p, B) + E T)/1 (p, B)] = 0, (65) 
with So(p,B) and £t, m (p, -B) given in Eqs. (3Z)-QS) and JHZJ-lEHl) respectively. The solution of JSBJ is 



In leading order in the Fermi coupling constant, expression l|66|) is approximated by 

E'i ~ |p| [l-a sin 2 a] - M.B±E (67) 

where 
and 

E = V2G F e-^ /2eB [7V - - JV+] (69) 

In the RHS of Eq. (|67|) , the first term is the modified neutrino kinetic energy, with the field-dependent radiative 
correction already found in the neutral case (see Eq. jQJl ); the second term can be interpreted as a magnetic- 
moment /magnetic-field interaction energy, with M playing the role of a neutrino effective magnetic moment; and the 
third term is a rest energy induced by the magnetized charged medium. 

We underline that, contrary to the case of the neutrino anomalous magnetic moment found beyond the standard 
model by including right-chiral neutrinos 30] , the induced magnetic moment here obtained does not require a massive 
neutrino. We remind that we are only considering left-chiral Dirac neutrinos, so the magnetic moment in l|67|) cannot 
be associated with the anomalous magnetic moment structural term a^F^ in the self-energy. This structure is 
forbidden in the present case by the explicit chirality of the standard model. In the case under study we have 
that although the neutrino is a neutral massless particle, the charged medium can endow it, through quantum 
corrections, with a magnetic moment proportional to the induced neutrino rest energy. Thus, we are in the presence 
of a peculiar magnetic moment created thanks to the particle-antiparticle unbalance of the charged medium. Such a 
charge asymmetry permits the formation of a net virtual current that, due to the magnetic field, circulates in a plane 
perpendicular to the neutrino propagation and therefore produces an effective magnetic moment in the direction of 
the neutrino momentum. 

In a charged medium CP-symmetry is violated |3l| , and a common property of electroweak media with CP violation 
is a separation between neutrino and anti-neutrino energies. In the present case this property is manifested in the 
double sign of Eq in (|67|l . Although this difference depends on the magnetic field entering in Nq, we stress that it 
depends on the magnitude of the magnetic field but not on its direction. Thus, the energy anisotropy connected to 
the first two terms of l|67|) is the same for neutrinos and anti-neutrinos. 

The neutrino/ antineutrino energy in the charged medium in the presence of a weak field has been found [15j to be 

E'J z = \p\-b'^±a', (70) 

with 

°' = ik {N - N ^ h ' = ik!wh^ {J -- M - (71) 

where N± are the electron/positron number densities, and f± the electron/positron distributions. 

Comparing the strong-field (|67(l and the weak-field i|70|) dispersion relations, we see that the field-dependent correc- 
tion to the kinetic energy appearing in the strong-field approximation is absent in the leading order of the weak-field 
case. For a strong-field, that correction turns out to be important for neutrinos propagating perpendicular to the 
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magnetic field. Another difference between (|67|l and J70|), is that in the last dispersion relation the rest energy a! 
does not depend on the magnetic field. This is the pure medium contribution that gives rise to the well known MSW 
effect 26]. As discussed above, in the strong-field case, where all parameters, including /; and |p| are assumed smaller 
than V eB, such a term is negligible compared to the contribution Eq which is proportional to the field. In both 
approximations there is an additional anisotropy related to the induced magnetic moment term, which is linear in 
the magnetic field and depends on the unbalance between particles and anti-particles. We should notice that this 
anisotropic term changes its sign when the neutrino reverses its motion. This property is crucial for a possible ex- 
planation of the peculiar high pulsar velocities [33 . Nevertheless, for the anisotropic term that modifies the kinetic 
energy in the strong-field approximation, the neutrino energy-momentum relation is invariant under the change of a 
by -a. 

Substituting (|67fl into (|59(l we obtain that the neutrino/anti-neutrino index of refraction in the charged medium is 
given by 

. 2 M-B T E 

n± = 1 + a sin a H — (72) 

The result Ij72(l implies that the index of refraction depends on the neutrino direction of motion relative to the magnetic 
field. Moreover, neutrinos and antincutrinos have different index of refractions and therefore different phase velocities 
even if they move in the same direction. For neutrinos, if Aq~ > (i.e. if there is a larger number of electrons than 
positrons) the index of refraction for B ^ is smaller than one, so their phase velocities will be smaller than light 
velocity. Thus, in the charged medium with strong magnetic field, neutrinos behave as massive particles. This is in 
agreement with the behavior they have in a dense medium, even in the absence of magnetic field |26j. On the other 
hand, for anti-neutrinos the index of refraction (the one with the positive sign in front of Eq in (J75J) can be larger 
than one. For instance, for anti-neutrinos moving opposite to the field-line directions the index of refraction is 

n_ = 1 + ^° > i (73) 

IpI 

Thus, antincutrinos will have phase velocities larger than c and depending on the magnetic field strength and 
electron density. In this regard we should mention that particles with zero rest mass can even have group velocities 
that exceed c in anomalously dispersive media |3^]. Moreover, in non- trivial backgrounds, particles with superluminal 
propagations are not intrinsically forbidden in quantum field theories. For instance, superluminal photons appears in 
curved spaces |34j. Casimir vacua [35| and in QED with compactified spatial dimensions |3r| . Discussions about the 
non- violation of microcausality by the existence of such superluminal velocities, as well as the lack of a contradiction 
between such a phenomenon with the bases of special theory of relativity, can be found in Refs. [23] , , • There, 
it was noted that the "front velocity," the one related to the index of refraction in the infinite frequency limit, is the 
one that cannot exceed c, since it is related to the signal transmission. 



VI. STRONG FIELD AND LLL APPROXIMATION. NUMERICAL TEST 

In the calculation of the neutrino self-energy in the strong-field limit (M\y VB ^> m e , T, /j,, |p|) done in Sec. IIVI 
we assumed that the main contribution to the self-energy loop diagram came from electrons in the LLL, since in the 
strong-field limit the energy gap between electron Landau levels is much larger than the electron average energy in 
the medium. In this Section we will check the validity of these arguments with the help of numerical calculations. 

Taking into account that the relevant physical quantity in this study is the neutrino energy, we can concentrate 
our analysis on the self-energy coefficients that give the largest contribution to the dispersion relation. It is easy to 
check, from the analytical structure of Eqs. 127|l - (|30|l . that the leading term in the expansion of the coefficients in 
powers of can be at most of order Keeping in mind that the only coefficients that could contribute 

with a 1/M^-order term in the dispersion relation JBJl are a and b (coefficients c and d appear squared in © and 
their contribution to the dispersion relation should start at least with terms of order 1/M^ V ) 1 we can approximate the 
dispersion relation © by 

^|p| 2 _ 2 ( a -&K (74) 
Then, to validate the LLL approximation we should numerically study the ratio 
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where a and b are respectively obtained from Eqs. 127|) and summing in all the Landau levels and taking 

parameter values in the strong- field region Mw 3> VB 3> m e: T, |p|. The coefficients ao and bo are the corresponding 
values at T = in the LLL . In the denominator of (|74() we are neglecting the thermal contributions ao(T) and bo (T) 
since they would make an insignificant contribution (1/Af^-order smaller) as compared with the vacuum ones in the 
parameter range T <C Mw- In Sec. IIVI we showed that b = 0, and that a is given by Eq. (|38|) . 

The strong-field conditions can be naturally found in many astrophysical systems like magnetars, neutron stars, etc. 
In cosmological applications, however, the viable primordial fields can never be much larger than the temperature 2 , 
as according to the equipartition principle, the magnetic energy can only be a small fraction of the universe energy 
density. This argument leads to the relation eB/T 2 ~ O(i). Clearly, this is not a situation very consistent with a 
strong-field approximation. However, even under this condition, it is natural to expect that the thermal energy should 
be barely enough to induce the occupation of just a few of the lower electron Landau levels, and one would not be 
surprised if the LLL approximation still gives the leading contribution to the dispersion relation. 

To investigate l|75[) it will be more convenient to independently study the variation range for each coefficient ratio, 

a(T = 0)+a(T) , , b(T = 0) + b(T) 

x a = - K —!--x a {B = Q,T = Q), K h = ± 1 K -^--xb(B = 0,T=0), (76) 

ao a 

where a(T = 0) (b(T — 0)) and a(T) (b(T)) are respectively the vacuum and thermal contributions of each coefficient 
a (6). Here we subtract the ultraviolet divergent zero-field, zero-temperature parts x a {B = 0,T = 0), x\,{B = 
0, T = 0), as they can only contribute, after renormalization, with negligible zero-field terms. The validity of the LLL 
approximation should imply that x a w 1 and y% -C 1 for the range of parameters considered. In the first part of this 
Section we examine the validity of the LLL in the zero-temperature case obtaining x a (T = 0) = a(T = 0)/ao — x a (B = 
0, T — 0) « 1, x;,(T = 0) = b(T = 0) /ao — Xb(B = 0, T = 0) <C 1. In the second part, we do a similar analysis, for the 
thermal part of a, finding that x a (T) = a{T)/ao <C 1, as expected from the theoretical considerations given in Sec. IIVI 
(see the discussion below Eq. I|48|l). As this last numerical calculation is performed for magnetic fields that can even 
be a few orders smaller than T 2 , the result x a (T) <gC 1 confirms the appropriateness of the LLL- approximation for a 
rather wide range of primordial fields. The analysis of the finite temperature part of coefficient b (xi>(T) — b(T)/ao) is 
not explicitly done in the paper, however, it is not hard to see that it gives rise to a similar result, that is, Xb(T) <C 1. 

A. Vacuum Contribution 

If we introduce in expressions l|27|l and (|28|l the integral representations 

1 



q + mj 



k 2 + M 2 



J da exp — (g 2 + m 2 )a, (77) 



d/3exp-(A: + M w )(3, (78) 



'vi- 
and take into account the recurrence relation |39j | 

(20 1 / a W-i(O = (^ + Ow(O, (79) 
in l|28|) . it is possible to rewrite coefficients a' and b' as 



. -iq 2 eB 

a = — R — 

2{2*f V \ 



<£x \ d A y / d 3 q / d 3 fc e -i(P^-Pi-Vx) e ^-(k+q-p)+iy-($'-k-q) 



J da J d/3{ e -^+™e)«-(fe!+^-eB )/3[(% _ q3){p3 +Pq) _ {qo + q3 )(p 3 - Po )r 2 t}S a Sp} (80) 





2 In cosmology the electric chemical potential is so small that we can always assume it is zero. 
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oo oo 

da ( d/3{e-^l^-(H+ M2 w-^[E+l p2 + l P_J^ pi]rtSa S f3 } (81) 
J 1 ~\~ ~t 1 £ 



Here, S a and Sp are found doing the sum in Landau levels with the help of Mehler's formula [38| 

S a = J2<Pl(P)Mp')t l = k(l - t 2 )]^ 1 / 2 exp - (^-) cxp - ^ ~ g (82) 



d 4 x / d 4 y / d 3 q / d 3 fc e -*CPi*i-pi-i/i) e *2-^+?-?)+ii?-^'-fc-?) 



where t — exp— (2eBa), p = \/2eB(xi + q^/eB) and p' — \ / 2eB(yi + q^/eB). Sp is obtained from Eq. H82|) by 
replacing t — ► r = cxp — (2eB(3), p — ► £ = V2e_B(xi — kz/eB) and — > £' = \[2eB{y\ — k 2 /eB) . 
Integrating in y, g and fc2, and introducing the variable changes xi — > £ and j/i — > we find 



oo oo oo oo oo 

, _ -ig 2 eBir 4 



„ = — ^ <5 4 (P-P') y d 2 k {l J daj d/3 J dt J de / { e -K«H-i'ii) a -H»a«-(*S+«^-^)/' e -*tt-C)pi/VS 

— oo — oo — oo 

[(1 + r 2 t)(p 2 -p 2 ) - (fco - fc 3 )(p 3 +Po) + (fco + fe 3 )(P3 -po^tJS^/j} (83) 

OO OO OO OO OO 

, = 2iff 2 (e£) 3 a 7T ^ (p _ j/) | d 2 h J da J dp J d ,J ^ {e - [ ( 9|| - P|| ) 2 + ^]a-(^ + A4-eB) /3 

— OO — OO — OO 

^ + + i^4 Pl + J- Pl J e-^-t'W^rtS'M (84) 



L l + i 1-t^ 1 V^B(l + <) 

where 



5 , = _ t2)] - 1/ 2 cxp _ ( K + p 2 /v^] 2 - £ ± w^] 2 } 



oyn [(r+P2/^)-(e+P 2 /^] 2 ^ 
cxp ^ 

After Wick rotating to Euclidean space and doing the Gaussian integrals in ka, ki, £ and we finally obtain 
*4/ ,s„2 2 1 ;a f 3 cosh(a + 2/3) 2 r>2 



o o 



cxp - [^^Fy + ,,, s ]_ Mh(f K ?±] (86) 
a + p 11 ctn(a) + cth(fj) 



oo oo 

6' = -<5 4 (p- P V^ 2 /** / d/3 r 4^-^S^exp(-m 2 a-M 2 ,/3) 

(a + p) smh (a + pj 
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-log m e /VeI 



FIG. 1: Plot of K a (T — 0) vs. logm e 1 for field range 10 2 m e < \feB < 10 2 m e and parameter values = |px| = 10 1 y / 2f, 
M w = 10 5 m e 




-log m e /VeI 



FIG. 2: Plot of X(,(T = 0) vs. logm e 1 for field range 10 2 m e < V eB < 10 2 m e and parameter values = \p±\ = 10 1 v2fi 
M w = 10 5 m e 



where we introduced the normalized parameters m e = m e /\/ eB, p^ = p^/y/eB, and Mw = Mw/VeB. Eqs. H86(l 
and (|87|l are as far as we can go in the calculation of a' and b' after summing in all Landau levels without using any 
approximation. 

In the zero-field limit the coefficients H86|) and H87J) reduce to 



oo oo 



a' B =o = b B=a = -^(P - P')9 2 ^ 2 J da J d/3 







3 , a 3 _ 7 —i . 

— — 3 exp-(a + — — -p M + M W 3) 
(a + 3) a + P 



(88) 



where Mw = Mw/m e , p^ = p M /m e . 

Our goal now is to investigate the validity of the LLL approximation for the zero-temperature contribution of the 
coefficients in the parameters' range: Mw ^ vl 3> m e , |p|. With this purpose, in Figs. 1 and 2 we plot, 



x a (T = 0) 



a - clb=o 
a 



vs logm" 1 G [-2,2] 



and 



K b (T = 0) = 6 vs logm- 1 G [-2, 2] 



no 



(89) 



(90) 
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respectively. Here ao is the LLL coefficient defined in Eq. I|38l) . Notice that primed and unprimed variables are related 
through: a' = (2tt) 4 5^ (p - p') a, b' = (2tt) 4 5^ (p - p') b. 

From Fig. 1 it can be seen that x a {T = 0) ~ 1 when y/eB > m e , |p| , while, from Fig. 2, we see that >r&(T = 0) <C 1 
in all the range 10 _2 rn e < V eB < 10 2 m e . This result corroborates the initial assumption that the LLL approximation 
gives the leading contribution to the vacuum part of the neutrino self-energy at strong fields (Mw ^> y/~B ^> m e , \p\). 



B. Thermal Contribution 



Let us consider now the finite temperature contributions. The finite temperature part of the coefficient a' can 
be found in the following way. We start from Eq. (JHHJl, perform the Wick rotation to Euclidean space, replace the 

oo 

integral in the four momentum by the temperature sum J dk^ — > (2nT) > (with k^ = 2irnT), and integrate in 

n— — oo 

&3 , to obtain 



, g 2 eBir^fTT 



oo oo oo oo 



a' = 8\p-p') y ^;' v " I da I dp I I dti'{S' a Sp 



— oo — oo 



y/a + P 



* — ' a + a 

n— — oo 



+*fc 4 [(P3 - Po)r 2 t - (p 3 + po )]} e -^+0)kl-2ia Po k i} 



After integrating in £ and and introducing the elliptic theta functions representation [36 

oo 

#3 (u/t) = exp[j7r(rn 2 + 2rra)] 

we find 



where 



n— — oo 



a ' = a (D +a (2) 



l« = <5 4 (p-p') 



(91) 



(92) 



OO oo 



x I da I dp sinhr 1 (a + P) ^ ^±^ 



exp - [p 2 + 77i 2 ]a exp — exp — 



Pi 







coth a + coth /? 



aP3 
a + P 



[r-H-^ipo+ps) + rt^(p 3 - Po )} + (pi - pfXr" 1 *- 1 /' + rt^) 



and 



3 (-2p Ta I i^f 2 {a + p) 



=5\p-p') 



(93) 



X dj dAsinh^A) 



i , , . eX ?iX^3 exp _ g| + falYt exp -M^(A - 7) exp - 



coth7 + coth(A — 7) 



d^ 3 (-2p T>y/i4nT*\ 

X [« 1/2 (p 3 -Pb) " ^H- 1/2 (pb +p 3 )] — 



In Eq. we introduced the notation 

r — exp — 2(A — 7), £ = exp — 2j 
Integrating (|94(l by parts and adding the result with l|93(l . we obtain 

/■ , .sinh-^A). 



a' = £>-p') : 



PoPn 



d\ — — [p 3 sinh(A) + p cosh(A)] [exp -[-pi + m 2 ]A] 6 3 [~2p T\ U\-kT 2 \ 



5 2 7r 2 0Fr /• rfA sinh- (A) [fe ginh(2A) + ^ cosh(2A)] r _m^ x 



PoPll 



A 



(0/i47rf 2 A) 



PoPu J J Va + j3 
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sinh 1 (a + /3) _ ^ 2 a ^ r _ 2 i _ 2l ^ or , 



■ exp— M w f3exp — [pu +m e ]aexp( -)p 



x exp - 



Pi 



coth a + coth , 



2ap 3 p 



- < ' [p 3 cosh(a + 2/3) + po sinhfa + 2/3)] - pfi cosh(a + 2/?) 

p y a + (3 " 



a + [3 sinh a sinh /3 [coth a + coth 0]' 
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To isolate the temperature part of a' , we have to subtract the zero-temperature (vacuum) contribution from 
(|96[) . With this aim, we use the Jacobi imaginary transformation [3^ 



(w/r) = (y/—ir) 1 exp[-inu 2 /T}e 3 ^- j 



to write 



6> 3 (-2p fa /i47rf 2 (a + /3) 



exp-[pgq 2 /(a + /3)] 
2T^^TW) 



1 + ' exp - 



n poan 



n— — 00 



4T 2 (a + /3) T(a + /3) 



where the symbol means that the term n = was taken out. 

n— — 00 

Then, to subtract the vacuum term from a' is equivalent to make in Eq. I|9tj|) the following substitution 

cxp-[p 2 a 2 /(a + /3)] 



• Ta /i47rT 2 (a + f3) } -> 9 3 ( -2p Ta i^T 2 (a + (3) - 



2T0V(a + (3) 
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log (a(T)/a ) 




FIG. 3: Plot of log(a(T)/oo) vs. log(f) for field range 1CT 2 T < y/eB < 10 2 T and parameter values |py | = \pj_\ = W^V^f, 
M w = 10 5 m e , m e = 0.1 
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(100) 
(101) 



(102) 



Eq. (|102fl gives the thermal part of the a' taking into account all Landau levels. To compare the thermal and the 
LLL contributions we plot loga(T)/ao vs. logT" 1 in Fig. 3. 
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From Fig. 3 we see that in all the range of temperatures considered, the thermal contribution is consistently smaller 
in several orders than the LLL contribution. This result, along with those obtained in Figs. 1 and 2, imply that 
thereby validating the LLL approximation in the strong and weakly-strong field cases {Mw ^> \ /f B ^> m e , |p|, 
eB > T 2 ). In addition, we would like to underline the consistency of these numerical results with the arguments used 
in Sec. IIVI (see discussion after Eq. I@5J)), in the sense that the finite temperature part should be much smaller than 
the vacuum contribution. 



VII. APPLICATIONS TO COSMOLOGY AND ASTROPHYSICS 

Two natural environments where the results of the present paper can be of interest are, on one side, stars possessing 
large magnetic fields (neutron stars, magnetars, etc), and, on the other, the early universe, which presumably was 
permeated by large primordial magnetic fields that eventually became the seeds , |4jj of the galactic fields that are 
observed today at scales of 100 Kpc J6J and larger. 

The effects of primordial magnetic fields on neutrino propagation are of special interest during the neutrino decou- 
pling era, on which a large number of these particles escaped the electroweak plasma with the possibility to develop 
significant oscillations. The strength of the primordial magnetic field in the neutrino decoupling era can be estimated 
from the constraints derived from the successful nucleosynthesis prediction of primordial 4 He abundance 0, as weu 
as on the neutrino mass and oscillation limits |4^ . These constraints, together with the energy equipartition principle, 
lead to the relations 

m 2 e < eB < ml, B/T 2 ~ 2, (103) 

with rrifi being the muon mass. Then, it is reasonable to assume that between the QCD phase transition epoch and 
the end of nucleosynthesis a primordial magnetic field in the above range could have been present 17]. 

Unlike the stellar material, whose density can be of order one, we know that the early universe was almost charge 
symmetric, with a particle-antiparticle asymmetry of only 10 -9 — 10 -10 . Therefore, when considering possible conse- 
quences of our results for the early universe, we should restrict the discussion to the neutral case (/i = 0). 

For the neutral case, it is known that at weak field all the corrections in the neutrino ener gy d ensity, whether they 
depend on the field and/or temperature, are second order in the Fermi coupling constant jl6|. [T^. [Til] , therefore 
negligible small. Nevertheless, it can be seen from Eqs. I|tjll) - (|62[) that if sufficiently strong magnetic fields were 
present in the primeval plasma, they would yield corrections to the energy density that are linear in the Fermi 
coupling constant. Moreover, as shown in Sec. IVII the strong-field approximation that led to these corrections is 
reliable even for magnetic fields in a more realistic range (i.e. as those satisfying condition l|103|) '). 

Notice that a field satisfying i|103|) would be effectively strong with respect to the electron-neutrino, but weak 
for the remaining neutrino flavors. If such a field existed during the decoupling era, it could significantly affect the 
Ve <-* Vp. , v T and v e <-» v s resonant oscillations 1 1 1 1] , as the field would differently modify the energy of v e compared 
to those of r/ft, v T and v s . The interesting new thing here is that despite these oscillations would take place in an 
essentially neutral medium, because of the strong field they will be as significant as those produced by the MSW 
mechanism [2^ in a dense medium. A peculiarity of the strong-field effects on the neutrino energy density is to give 
rise to anisotropic resonant oscillations. That is, the oscillation probability depends on the direction of the neutrino 
propagation with respect to the magnetic field. 

Another interesting question related to primordial magnetic fields is whether they influenced neutrino propagation 
prior to the electroweak phase transition, since some of the mechanisms of primordial magnetic field generation allow 
their existence at very early epochs 0, @. Before the electroweak transition a primordial magnetic field could 
only exist in the form of a U(l) hypermagnetic field 45]. Any non-Abelian "magnetic" field would decay at high 
temperatures because it would acquire a non-perturbative infrared magnetic mass g 2 T. The implications of primordial 
hyp ermagnetic fields in neutrino propagation before the electroweak phase transition have been studied in Refs. |43| . 
|44j . In [43 the strong-field effects led to a large anisotropy in the leptons' energies. The anisotropy was due to the 
degeneracy in the energy, which in the leading order does not depend on the transverse momentum component. It is 
in contrast with the magnetic field case, which takes place after the electroweak phase transition, where, although an 
anisotropy in the neutrino energy is also present, no degeneracy occurs (see Eas. I61|l - H62f) '). This different behavior 
is a direct consequence, on one hand, of the nonzero neutrino hypercharge and hence of the minimally coupling of 
these particles to the hypermagnetic field; and on the other, of their electrical neutrality, which implies that they can 
interact with a magnetic field through radiative corrections only. 

We underline that the anisotropy in the energy spectrum of neutrinos in background fields could provide an inde- 
pendent way to verify the existence of primordial magnetic fields, since a field-induced anisotropy would be reflected 
as an imprint in a yet undetected and elusive cosmic background of neutrinos produced during the decoupling era. 
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For astrophysical applications we should turn our attention to the star interiors characterized by high densities 
of charged particles under strong magnetic fields. As shown in Sec. in the charged medium case the strong-field 
approximation gives rise to a modification of the kinetic part of the neutrino energy (first term in Eq. (J^TJ ) , that is not 
present in leading order at weak field (Eq. J7UJ0 15]. Due to this new term the largest field-dependent contribution 
to the neutrino kinetic energy comes from neutrinos propagating perpendicularly to the field. 

It is worth noticing that the anisotropic term associated with the effective magnetic moment in the dispersion 
relations (|67l) and l|7U|) changes sign when neutrinos reverse their direction of motion. On the other hand, the kinetic 
energy term in (|67H remains unchanged when the neutrinos moves in opposite direction, even though it depends on the 
direction of propagation of the neutrinos with respect to the magnetic field. As known [32} , the anisotropy associated 
to the magnetic-moment contribution to the neutrino energy can be relevant for a possible explanation of the peculiar 
high pulsar velocities. It would be interesting to consider the combined effects of the two different anisotropics to 
understand whether they could affect the dynamics of proto-neutron stars. 

Another possible ground of applications of our results is supernova neutrinos. Core collapse supernovae are dom- 
inated by neutrinos and their transport properties. In addition, the observation of neutrinos from supernova, which 
is essential to confirm the basic picture of supernova explosion, can be affected by neutrino oscillations. As we have 
pointed out in this paper, neutrino transport and oscillations can be both modified by the presence of a strong mag- 
netic field. Magnetic fields as strong as 10 14 to 10 16 G could exist in the first seconds of neutrino emission inside the 
supernova core |46j |. Thus, the electron- neutrino energy spectrum found in this paper, within the strong- field limit 
for the charged medium, should be considered for any study of neutrino oscillations under those conditions. 

VIII. CONCLUDING REMARKS 

In this paper we carried out a thorough study of the propagation of neutrinos in strongly magnetized neutral and 
charged media. We started from the most general structure of the neutrino self-energy in a magnetic field, expressing 
it as the sum of four independent covariant terms with coefficients that are functions of the physical variables of the 
theory and whose values depend on the approximation considered. General expressions of the four coefficients at 
one-loop approximation were given in Eqs. (|27I30|) . 

The coefficients were then calculated in the strong-field limit using the LLL aproximation for the electrons. The 
LLL was assumed to be valid in the parameter range M w > V~B > m e , |p|, eB > T 2 . To justify it one should keep 
in mind that under these conditions most electrons would not have enough energy to overcome the gap between the 
Landau levels. Hence, they will be mainly confined to their lower levels and the leading contribution would come from 
the LLL. This assumption was also corroborated for the above parameters' range by numerical calculations summing 
in all Landau levels. 

The dispersion relation of the neutrinos was written as a function of the four coefficients of the self-energy structures, 
allowing in this way to straightforwardly obtain the neutrino's energy in the strong-field limit for each physical case. 

In concordance with results previously obtained in charged media at weak fields [l5j|,^4|], m the strong- field case an 
energy term associated with the interaction between the magnetic field and the effective magnetic moment was also 
found at leading order in Gf- This interaction energy disappears in the neutral medium, since in a charged-symmetric 
plasma the contribution to the effective magnetic moment coming from electrons and positrons cancels out. 

A main outcome of our investigation was to show that in strongly magnetized systems a term of different nature 
emerges in both charged and neutral media. The new term, which is linear in the magnetic field and of first order in 
Gf, enters as a correction to the neutrino kinetic energy in the presence of a strong-magnetic field. This correction is 
present even in a strongly magnetized vacuum, since it is related to the vacuum part (T = 0, fi = 0) of the neutrino 
self-energy at B ^ 0. 

A characteristic of the field-dependent corrections to the neutrino energy is that they produce an anisotropic index 
of refraction, since neutrinos moving along different directions have different field-dependent dispersion relations. We 
should underline that while the magnetic moment interaction term produces a maximum field effect for neutrinos 
propagating along the field lines, the field correction to the kinetic energy does not contribute to those propagation 
modes, but on the contrary, the maximum kinetic-energy effect takes place for neutrinos propagating perpendicularly 
to the field direction. We stress that the anisotropy does not differentiate between neutrinos and antineutrinos. 

The charged medium results reported in the current work can be of interest for the astrophysics of neutrinos in stars 
with large magnetic fields. On the other hand, our finding for the neutral medium can have applications in cosmology, 
if the existence of high primordial magnetic fields is finally confirmed. Contrary to some authors' belief that, 
regardless of the field intensity, the neutrino dispersion relation in the early universe is well approximated by the 
dispersion relation in the zero field medium, our results indicate that strong, and even weakly strong, magnetic fields 
can give rise to a contribution to the neutrino energy that is several orders larger than the pure thermal contribution. 

The field-dependent correction to the neutrino energy in a neutral medium with strong magnetic field can have 
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an impact in neutrino flavor-oscillations in the primeval plasma |ll| and therefore affect primordial nucleosynthesis. 
Hence, this new effect could be important to establish possible limits to the strength of the primordial magnetic field. 
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